Phenomenological Ginzburg-Landau-like theory for superconductivity in the cuprates 
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We propose and develop here a phenomenological Ginzburg-Landau-like theory of cuprate high- 
temperature superconductivity. The free energy of a cuprate superconductor is expressed as a 
functional T of the complex spin-singlet pair amplitude ipij = ip m = A m exp(i^> m ) where i and 
j are nearest-neighbor sites of the square planar Cu lattice in which the superconductivity is be- 
lieved to primarily reside and m labels the site located at the center of the bond between i and j. 
The system is modeled as a weakly coupled stack of such planes. We hypothesize a simple form, 



•HA, <t>] = E m {AM> + (B/2) A 4 m ) + C £ <mn> A m A n cos( 



), for the functional, where m and 



n are nearest-neighbor sites on the bond-center lattice. This form is analogous to the original con- 
tinuum Ginzburg-Landau free energy functional; the coefficients A, B and C are determined from 
comparison with experiments. A combination of analytic approximations, numerical minimization 
and Monte Carlo simulations is used to work out a number of consequences of the proposed func- 
tional for specific choices of A, B and C as functions of hole density x and temperature T. There 
can be a rapid crossover of (A m ) from small to large values as A changes sign from positive to 
negative on lowering T; this crossover temperatures T lns (x) is identified with the observed pseudo- 
gap temperature T*(x). The thermodynamic superconducting phase-coherence transition occurs at 
a lower temperature T c (x), and describes superconductivity with d-wave symmetry for positive C. 
The calculated T c (x) curve has the observed parabolic shape. The results for the superfluid density 
p s (x,T), the local gap magnitude (A m ), the specific heat C v (x,T) (with and without a magnetic 
field) as well as vortex properties, all obtained using the proposed functional, are compared success- 
fully with experiments. We also obtain the electron spectral density as influenced by the coupling 
between the electrons and the correlation function of the pair amplitude calculated from the func- 
tional and compare the results successfully with the electronic spectrum measured through Angle 
Resolved Photoemission Spectroscopy (ARPES). For the specific heat, vortex structure and elec- 
tron spectral density, only some of the final results are reported here; the details are presented in 
subsequent papers. 



I. INTRODUCTION 



The last two decades have seen unprecedented ex- 
perimental and theoretical activities involving cuprates 
which exhibit high-temperature superconductivity [1-4]. 
Even after this long period of research which has seen 
dramatic advances in experimental techniques and mate- 
rials quality, as well as discovery of many unusual phe- 
nomena such as the ubiquitous pseudogap in undcrdoped 
cuprates [5-8] and the 'strange metal' phase above the 
superconducting transition temperature around optimal 
doping [1-3], there is no common, broadly accepted un- 
derstanding yet about their origin. 

Motivated by the above, especially the increasing vol- 
ume of sophisticated spectroscopic data on the cuprates 
(such as those obtained from ARPES [9, 10], STM [11] 
and Raman [12] experiments), we propose and develop 
here, as well as in subsequent papers, a new phenomeno- 
logical model for cuprate superconductivity that is anal- 
ogous in form to the well-known Ginzburg-Landau (GL) 
theory [13] of superconductivity. The starting point of 
our description is the assumption that the free energy 
of a cuprate superconductor can be expressed as a func- 
tional solely of the complex pair amplitude. In the orig- 
inal continuum GL theory, the free energy, expressed 
as a functional of the complex order parameter field 



ip{r) = A(r) exp (i</>(r)), has the form 
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This form is justified near the actual superconducting 
transition where the magnitude of the order parameter is 
small, so that a low-order power series expansion in ?/>(r) 
is adequate. Further, ip(r) is assumed to vary slowly with 
r so that it suffices to keep only the | V-0(r)| 2 term; this is 
the case in conventional superconductors where the nat- 
ural superconducting length scale (also the coarse grain- 
ing scale) £o is large (compared, say, to the Fermi wave- 
length). After the advent of the microscopic Bardeen- 
Coopcr-Schricffcr (BCS) [14] theory of superconductivity, 
-0(r) was identified by Gor'kov [15] with the Cooper pair 
amplitude, i.e. V( r ) = ( a t( r ) a 4-( r ))' where- a a {r) (aj^r)) 
is the operator which destroys (creates) an electron at r 
with spin a (a =t,J,). Gor'kov also obtained the coeffi- 
cients A c , B Cl C' c in terms of the electronic parameters 
of the metal. 

In our phenomenological description, we hypothesize 
that a free energy functional similar in structure to that 
of Eq.(l), but defined on the square planar Cu02 lattice, 
describes the properties of cuprate superconductors for a 
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FIG. 1: The square Cu lattice sites i,j,k,l,.. in the CuCb 
plane and construction of the bond lattice out of the centers 
of the Cu-O-Cu bonds. The solid circles at {Hi = i} (blue) 
represent the positions of Cu lattice sites and {R m = m = ij} 
(magenta) the positions of bond centre lattice sites. Alterna- 
tively, we denote the bond centre lattice site between Ri and 
Rj — Hi + aft as R; M = Hi + (a/2)fi with ft — +x, +y. The 
arrows indicate the direction of equivalent planar spins, with 
S m = (A m cos 4> m , A m sin <j>m) representing the complex order 
parameter ipij = ip m = A m exp(i</f) m ) and antiferromagnetic 
ordering (shown) of spins translating into a d-wave symmetry 
gap (long-range order). 



FIG. 2: A schematic illustration of the hole doping x and 
temperature T plane (entire shaded region) where we assume 
the functional of Eq.(2) to be applicable. Tf(x) (solid brown 
line) and T c (x) (solid blue line) are shown along with the 
experimental superconducting (SC) dome and antiferromag- 
netic (AF) regime at very low hole doping. The two arcs 
shown by dotted lines denote regions where quantum fluctu- 
ation effects, as well as other low-energy degrees of freedom, 
such as electronic and spin plus their coupling with pair de- 
grees of freedom, need to be explicitly included in the free 
energy functional. For instance, inclusion of quantum phase 
fluctuation effects in a minimal level leads to a T c (x) curve in 
agreement with experiment (See Section III). 



fairly wide range of hole doping (x) and temperature (T). 
Fig.l shows the square planar lattice schematically, and 
Fig. 2 the region of the (x, T) plane where our phenomeno- 
logical description is assumed to be applicable. The free 
energy is assumed to be a functional of the complex spin- 
singlet pair amplitude ipy = ip m = A m exp(i<p m ) where i 
and j are nearest-neighbor sites of the square planar Cu 
lattice and m labels the 'bond-center lattice' site located 
at the center of the bond between the lattice sites i and j 
(see Fig. 1). The highly anisotropic cuprate superconduc- 
tor is modeled as a weakly coupled stack of Cu02 planes 
in which the superconductivity is believed to primarily 
reside and wc ignore, as a first approximation, the inter- 
plane coupling. The free energy functional for a single 
Cu02 plane is assumed to have the form 

F{{ A m , <t> m }) = T ({A m }) + Ji({A m , 4> m }), (2a) 

Jo({A„J) = J2 ( A ^n + f Al) , (2b) 

Jl({A m , 4> m }) = C ^ A m A„ COs(0 m - 0n)(2c) 



A Gor'kov like interpretation of tpij is that it is the aver- 
age spin-singlet nearest-neighbor Cooper pair amplitude, 
i.e. i/iij = (bij)/V2 = (l/2)(ai4.ajt - a^aji). The sites 
i and j are different because strong electron repulsion 
(symbolized for example by the Mott-Hubbard U) disfa- 
vors on-site pairing, while the existence of large nearest- 



neighbor antiferromagnetic spin-spin interaction in the 
parent cuprate is identically equivalent for spin-^ elec- 
trons to attraction between nearest-neighbor pairs (i.e. 
Jij(Si.Sj — ifiifij) = —Jijbjjbij with Si and hi the elec- 
tron spin and number operators respectively at the z-th 
site) . This favors the formation of nearest- neighbor spin- 
singlet pairs. 

The first part Tq of T is the sum of identical inde- 
pendent terms each of which is a function of only the 
magnitude A m of the order parameter on the bond lat- 
tice site. Eq.(2b) is a simple form for it in the image of 
Eq.(l), with A and B depending, in general, on x and 
T . We assume that B is a positive constant independent 
of x and T and choose A(x,T) to change sign along a 
straight line T®(x) running from T = Tq at x = to 
T = at x = x c (see Fig. 2). As a first approximation, 
this line can be identified with the pseudogap tempera- 
ture T* (x) because the magnitude of the local pair am- 
pliitude, (A m ), can increase dramatically as the temper- 
ature crosses this line, so that A changes from a positive 
to a negative value. The occurrence of superconductivity, 
characterized by a nonzero stiffness for long- wavelength 
phase fluctuations and the associated superconducting 
phase coherence, depends on the phase coupling term, 
Eq.(2c). If C in Eq.(2c) is taken to be proportional to x, 
the superconducting transition temperature T c , as calcu- 
lated in our theory, turns out to be proportional to x for 



small values of it, in conformity with what is observed, 
e.g. the Uemura correlation [16]. Also, if C is taken to be 
positive, the transition is to a d-wave symmetry super- 
conducting state (see Section II). We, therefore, make 
this choice. 

We emphasize that the assumed form of the func- 
tional and the dependence of the coefficients on x and 
T are purely phenomenological, guided by experimen- 
tal results - the functional is not derived from a micro- 
scopic theory. The functional satisfies the usual symme- 
try and stability requirements: the absence of odd pow- 
ers of ipm ensures invariance of the free energy under a 
global change of phase, and the free energy is bounded 
below for the chosen positive B. Since A m = \ip m \ and 
A m A„cos(0 m - </>„) = -Qip m - ij n \ 2 - A 2 m - A*)/2, 
it is readily seen that the free energy of Eq.(2) is simi- 
lar in form to a discretized version of the GL functional 
of Eq.(l). However, there are important differences be- 
tween our phenomenological approach and the original 
GL theory - these differences are discussed in detail in 
Section II. 

The main objective of our study is to investigate 
whether the free energy functional defined above pro- 
vides a good description of experimental results over a 
wide range of x and T. To this end, we have carried 
out several investigations of the thermodynamic behav- 
ior of a system whose equilibrium properties are given 
by canonical (thermal) averages with the functional of 
Eq.(2) playing the role of the 'Hamiltonian' or energy 
function. These calculations have been performed at sev- 
eral levels of sophistication. We first used simple single- 
site mean-field theory to obtain qualitative information 
about the behavior of the system over a wide range of x 
and T. We also used cluster mean-field theory to obtain 
more accurate estimates of the superconducting transi- 
tion temperature as a function of doping. We used nu- 
merical minimization of the free energy to obtain exact 
results for the properties of the system and the structure 
of vortices at zero temperature. We also used extensive 
Monte Carlo (MC) simulations to obtain exact (modulo 
finite-size effects) information about the thermodynamic 
behavior of the system at finite temperatures. Since the 
free energy of Eq.(2) may be viewed as the Hamiltonian 
of a two-dimensional XY model with fluctuations in the 
magnitudes of the 'spins' (see Section II for the details of 
this analogy) , we made use of well-known results about 
the behavior of the XY model in two dimensions in the 
analysis of the data obtained from our MC simulations. 
Finally, we extended our free energy functional to include 
quantum phase fluctuations (see Section III) in order to 
study the effects of these fluctuations on the transition 
temperature, and included coupling of the pair degrees 
of freedom to electrons (see Section VIII) to study the 
spectral properties of electrons measured in ARPES ex- 
periments. Simple, physically motivated, approximate 
analytic methods were used in these studies. The main 
results obtained from these extensive analytic and nu- 
merical calculations are summarized below. 



As a starting point, we calculate the superconducting 
transition temperature T c (x) and the average magnitude 
of the local pair amplitude, (A m ), using single-site mean- 
field theory for the model of Eq.(2). We show that this 
approximation leads to general features of the x~T phase 
diagram in agreement with experiment. In particular, we 
find a phase coherent superconducting state with d-wave 
symmetry below a parabolic T c (x) dome and a phase in- 
coherent state with a perceptible local gap that persists 
up to a temperature around T®(x). Further, effects of 
thermal fluctuations beyond the mean-field level are cap- 
tured via MC simulations of the model of Eq.(2) for a 
finite two-dimensional lattice. Section III describes the 
results for T c (x) obtained from these simulations. The 
actual values of A, B and C used in these calculations 
are discussed in Section II. 

The superfluid stiffness p a (x,T) (a quantity measured 
e.g. via the penetration depth) is calculated in Section 
IV. Its doping and temperature dependence compare well 
with experimental results [17-21]. The thermally aver- 
aged local gap A(x,T) = (A m ) is obtained in Section V 
where we calculate the temperature T ms (x) correspond- 
ing to the maximum slope of this quantity with and with- 
out the C term. This temperature provides a measure of 
the pseudogap temperature T*{x). We use these results 
to remark on contrasting scenarios [7, 8] proposed for the 
doping dependence of the pseudogap. We find that there 
is a contribution to A(x,T) that 'turns on' at T c (x), the 
superconducting transition temperature. This is obvi- 
ously connected with persistent observations of two dif- 
ferent kinds of energy gaps in several experiments [22, 23]. 
We also calculate the ratio 2A(x, 0)/T c (x) which is ob- 
served to be generally much larger than the BCS value 
of about 4 over a wide range of a; [11, 24], and to vary 
from system to system within the cuprate family for the 
same x. Our results rationalize this behavior, which is 
expected here since the origins of A(x, 0) and T c are dif- 
ferent. 

The contribution of the pair degrees of freedom to ther- 
mal properties, such as the specific heat C v , can be ob- 
tained from the free energy functional of Eq.(2). We 
briefly report in Section VI our calculation of C v (details 
are given in a subsequent paper [25]), and find that there 
are two peaks [26-29] in it, a sharp one connected with T c 
(ordering of the phase of %p m ) and a relatively broad one 
('hump') linked to T* (rapid growth of the magnitude of 
-0m)- The former is specially sensitive to the presence of 
a magnetic field, as we find in agreement with experiment 
[30, 31]. Vortices, which are topological singularities in 
phase, are naturally explored in our approach [32]. We 
use the functional of Eq.(2) to find A m and </> m at differ- 
ent sites m for a 2tt vortex whose core is at the center of 
a square plaquette of Cu lattice sites (Section VII). We 
find that the vortex changes character from being pri- 
marily a phase or Josephson vortex for small x to a more 
BCS-like vortex with a large diminution in the magni- 
tude A m as one approaches the vortex core for large x. 
Rcf. [33] describes these results in greater detail. 



Experimental information about the pair field i\) m and 
its correlations is not obtained directly, but from its cou- 
pling to electrons (e.g. ARPES [9, 10] and STM [If]), 
photons (e.g. Raman scattering [f 2] and light absorption 
[34] ) and neutrons [4] . We therefore develop a theory for 
the coupling of electrons near the Fermi energy with ip m 
and outline it in Section VIII. A separate paper [35] de- 
scribes this approach in detail as well as the results (e.g. 
Fermi arcs that are ubiquitous above T c , and the pseudo- 
gap for various momentum regions of the Fermi surface, 
especially the antinodal region) which compare very well 
with the results of recent ARPES measurements. We 
present here the results for the antinodal pseudogap fill- 
ing temperature T an (x) and compare it with the other 
estimate of the pseudogap temperature, T ms (ir), obtained 
in Section V. 

The results summarized above establish that the sim- 
ple free energy functional proposed here provides a con- 
sistent and qualitatively correct (quantitative in some 
cases) description of a variery of experimentally observed 
properties of cuprate superconductors over a wide range 
of temperature T and hole concentration x. This is the 
main conclusion of our study. Section IX discusses cer- 
tain generalizations, applications and limitations of the 
approach used in our study. Appendices A, B and C 
describe some technical details of the calculations. 



II. THE FREE ENERGY FUNCTIONAL 
A. Generalities 

As noted above, the free energy functional used in our 
study is phenomenological in nature with experimentally 
inspired coefficients. We have deliberately kept it as sim- 
ple as possible, without violating basic requirements of 
symmetry and stability. The form of the functional of 
Eq.(2) is analogous to that used in conventional GL the- 
ory. However, our approach is different from the GL 
theory in several ways. The form of the free energy func- 
tional used in the GL theory of supcrcoductivity and in 
similar theories of other continuous phase transitions [36] 
can be justified only if the temperature is close to the 
transition temperature. This approach, therefore, is ex- 
pected to yield quantitatively correct results only in the 
vicinity of the superconducting transition. This regime of 
validity is ordained by the requirement of smallness and 
slow spatial variation of the order parameter. Our use 
of the simple, GL-like functional of Eq.(2) over a broad 
(x, T) region can not be justified from similar consider- 
ations: the validity of our approach can only be judged 
a posteriori by comparing its consequences with exper- 
iments. Hence, we have calculated a variety of exper- 
imentally measurable quantities using the functional of 
Eq.(2) and compared the results with those of experi- 
ments. As discussed in detail in subsequent Sections, we 
find qualitative (and quantitative in some cases) agree- 
ment between the theoretical and experimental results for 



a wide variety of properties of cuprate superconductors. 
This establishes the usefulness of our phenomenological 
approach in describing the properties of cuprate super- 
conductors over a wide range of x and T. 

Another important difference between our approach 
and conventional GL theory is that the free energy func- 
tional we consider is not coarse-grained in the GL sense. 
We believe that this is natural because all cuprate su- 
perconductors are characterized by short intrinsic pairing 
length scales or coarse-graining lengths (£o ~ 15 — 20 A in 
the cuprates rather than the value of ~ 10, 000 A for 'con- 
ventional' pure superconductors) . We thus use a 'nearest- 
neighbor' coupling of the pair amplitudes defined at the 
sites of the atomic bond lattice in the second term of 
our functional (Eq.(2c)). Another difference between the 
functional used in our study and that of conventional 
GL theory is that the sign of the coupling constant C in 
Eq.(2c) is taken to be positive, so that the pair ampli- 
tudes at nearest-neighbor sites of the bond-center lattice 
have a phase difference of 7r in the ground state. This dif- 
ference in sign between the pair amplitudes on the 'hor- 
izontal' (in the x-dircction) and 'vertical' (y-direction) 
bonds of the Cu lattice corresponds to the superconduct- 
ing state having d-wave symmetry. This is consistent 
with the experimental fact that the superconducting gap 
Ak is proportional to (cos k x a — cos k y a), which arises in 
our description from a combination of nearest-neighbor 
Cooper pairs with relative phases as mentioned above. 

Some of the methods of calculation used in our study 
are also different from that in the conventional GL theory 
of superconductivity in which physical properties are cal- 
culated using simple mean field theory. The mean-field 
results are expected to be valid if the temperature is out- 
side the so-called 'critical' region [36] around the tran- 
sition temperature where the effects of fluctuations, not 
included in a mean-field analysis, are important. For con- 
ventional superconductors with long coherence lengths, 
the width of the critical region is very small, so that 
mean field theory provides a good description of most of 
the experimentally observed behavior. This, however, is 
not the case for cuprate superconductors with very short 
coherence lengths and for our model of cuprate supercon- 
ductivity. For this reason, we have to go beyond mean 
field theory (which provides a qualitatively correct, but 
not quantitative description of the general behavior) and 
use other methods (such as MC simulations) to obtain 
accurate results for the thermodynamic behavior of our 
model. 

A natural description of the pair amplitude ip m is as 
a planar spin of length A m pointing in a direction that 
makes an angle 4> m with a fixed axis. The thermal (Boltz- 
mann) probability of the length distribution is given pri- 
marily by J ({Am}) of Eq.(2b) and the term in Eq.(2c) 
may be thought of as the coupling between such 'spins'. 
The temperature T* (x) can be identified roughly as that 
at which the 'spin' at each bond lattice site acquires a 
sizable length locally without any global ordering of the 
angles, whereas the 'antiferromagnctic' (C > 0) nearest- 



neighbor interaction leads to global order (d-wave super- 
conductivity) setting in at T c . The two temperatures are 
well separated for small x because A, B and C are so 
chosen that T*(x ~ 0) >> T c . The region between T* 
and T c is the pseudogap regime where in the spin lan- 
guage, antiferromagnetic short-range correlations grow 
with decreasing temperature, its length scale diverging 
at T c . There is considerable experimental evidence for 
this view [5-7], though there is also the alternative view 
that T*(x) is associated with a new long-range order, e.g. 
d-density wave (DDW) [37] or time reversal symmetry 
breaking circulating currents [38]. 

The BCS theory and conventional GL theory in which 
the 'spin' formation and ordering temperatures are the 
same are limiting cases of this scenario. Something like 
this is expected to happen in cuprates near x c (Fig. 2) 
as also follows from our functional. The state below T c 
has nonzero order parameter (ipm) for a system above 
two dimensions, and is a Berezinskii-Kosterlitz-Thouless 
(BKT) [39-41] bound vortex state with quasi long range 
order in two dimensions, in which case T c is identified 
with the vortex unbinding temperature Tbkt- In the 
former case, the order parameter is the sublattice mag- 
netization Ad(x,T) = \(ipm}\ with a k-dependent gap 
Ak = (Ad/2)(cos k x a — cos k y a) . The interlayer coupling 
can be described, a la Lawrence and Doniach [42], by 
adding say a nearest-neighbor coupling between 'spins' 
on different layers to our functional in Eq.(2). Since this 
is in practice relatively small (the measured anisotropy 
ratio in Bi2212 is about 100, for example [43]), it makes 
very little difference quantitatively to most of our es- 
timates which generally neglect this coupling. For in- 
stance T c calculated by estimating the BKT transition 
temperature (Tbkt) from MC simulation of the two- 
dimensional model of Eq.(2) (See Section III) is expected 
to be very close to the actual transition temperature in 
the anisotropic 3D model with such small interlayer cou- 
pling. 

A conventional GL theory of cuprate superconductiv- 
ity would involve a functional similar to that in Eq.(l) 
(but with additional terms allowed by symmetry) with 
"0(r) the d-wave superconducting order parameter, and 
the coefficient so chosen that a mean-field treatment of 
the free energy leads to a dome-shaped T c (x) curve sim- 
ilar to that found in experiments. However, a mean-field 
treatment and the conclusions obtained from it would 
not be reliable because of the smallness of the supercon- 
ducting coherence length in the cuprates and consequent 
large fluctuation effects. In particular, the pseudogap 
temperature T* (which is much larger than T c for small 
x, and goes to T c as x increases) would be absent in 
such a theory. By contrast, we assume here that the 
basic low-energy Cooper pair degree of freedom in the 
cuprates is the bond pair, give a physical meaning to T* 
as a pair magnitude crossover temperature, and describe 
the regime between T* and T c as one in which the correla- 
tion length associated with superconducting fluctuations 
of d-wave symmetry grows and diverges at T c . The effect 



of these fluctuations is found to be crucial for many phys- 
ical properties, e.g. the Fermi arc phenomenon, and the 
filling of the antinodal pseudogap as T rises to T* . The 
superconducting order with d-wave symmetry that sets 
in at T c is an emergent collective effect, arising from the 
short-range ipm^n interaction, much as long-range Nccl 
order arises from an antifcrromagentic coupling between 
nearest-neighbor spins. 

GL theories for cuprates have been proposed by a 
large number of authors, arising either out of a particu- 
lar model for electronic behavior and often coupled with 
the assumption of a particular 'glue' for binding electrons 
into pairs [44-46] , or out of lattice symmetry considera- 
tions [47, 48]. The functional in Eq.(2) is consistent with 
square lattice symmetry and, in principle, does not as- 
sume any particular electronic approach (weak coupling 
or strong correlation, for example) or a mechanism for the 
'glue'. However, some of the properties of the coefficients 
are natural in a strong electron correlation framework. 
For example, mobile holes in such a system can cause 
a transition between a state in which there is a Cooper 
pair in the x directed ij bond (Fig.l) to one in which the 
Cooper pair is in an otherwise identical but y directed 
bond jk nearest to it (or vice versa), thus leading to a 
nonzero term T\ in Eq.(2). This is probably connected 
with the observed [49] empirical correlation between T c 
and the diagonal or next-nearest-neighbor hopping am- 
plitude of electrons in the Cu lattice. 



B. Parameters of the Functional 

The coefficients A, B and C are chosen to be consis- 
tent with experiments. Specifically, the coefficients are 
as follows: 



x,T) 
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T -Toll- 
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e T/T p ^ 


(3a) 


B 


= b t , 






(3b) 


C(x) 
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We require A m to have dimensions of energy [E] (or tem- 
perature for Boltzmann constant ks = 1) and hence Aq, 
Bo and Co have dimensions of [E]~ 2 , [E]~ 4 and [E}~ 2 
respectively. They are rewritten in terms of To as well 
as three dimensionless parameters /, b and c so that 
T carries dimension of energy as well. We thus have, 
Ao = (//To) 2 , B = b(f/T ) 4 and C = c(//T ) 2 . We 
choose b and c to have values of order unity and fix them 
for different hole doped cuprates by comparing Ao(x), 
T*(x) and T° pt obtained from the theory with experi- 
ments (see below for details). 

The two temperature dependent parts of A as given 
above arise as follows. The part [T — To(l — x/x c )] reflects 
our identification of the zero of A(x, T) with the pseudo- 
gap temperature and the experimental observation that 
the pseudogap region extends downwards nearly linearly 



from T = T at x = to T = for x = x c . The rela- 
tion between this straight line Tj (x), the experimental 
T*(x) and the related quantities T^(x) (obtained from 
a maximum slope criterion, Section V) as well as T an (x) 
(obtained from the antinodal gap filling criterion for the 
electron spectral function, Section VIII) is shown in Fig. 7 
and Fig. 16. The exponential factor e T / T p suppresses 
A(x,T) at high temperatures (T >> Tf{x)) with re- 
spect to its temperature independent cquipartition value 
y/T/A(x, T) which will result from the classical func- 
tional (Eq.(2)) being used well beyond the near prox- 
imity of any critical temperature where it is valid. Such 
a suppression is natural in a degenerate Fermi system; 
the relevant local electron pair susceptibility is rather 
small above the pair binding temperature and below the 
degeneracy temperature. The temperature scale T p is 
of order To, this being the energy scale for pair bind- 
ing. We take it to be To unless stated otherwise. In all 
the calculations below, we choose x c — 0.3 and b = 0.1 
(except in Fig. 4(b)). b along with T p controls the tem- 
perature dependence of A(x,T), especially the decrease 
of A(x, T) across the 'pseudogap temperature' line T*(x) 
and other details such as the height of the specific heat 
hump around T*(x). Values of /, c and To can be fixed 
for a variety of cuprates by comparing zero temperature 
gap Aq(x), T*(x) and T° pt with experiments. For exam- 
ple, a choice of parameters, roughly suitable for Bi2212, 
which has an experimental T° pt ~ 91 K, gives / ~ 1.33, 
c ~ 0.3 with A (x = 0) ~ 82 mcV, T ~ 400 K and 
T°p* t ~ 72 K (T c °p* ~ 110 K from single site mean field 
theory, see Section III). Unless otherwise stated, we have 
used the above choice of parameter values in the rest of 
the paper. 



III. SUPERCONDUCTING TRANSITION 
TEMPERATURE T c (x) 

The superconducting state is characterized by macro- 
scopic phase coherence. For superconductivity in 
cuprates described by the functional (Eq.(2)) this means 
a nonzero value for the superfluid stiffness or supcrfluid 
density p s (x,T) given by the formula [50], 
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(4) 



where the subscript m + \i refers to R m + Ip, with p run- 
ning over x and y directions in the bond lattice coordinate 
system (rotated by 45° with respect to the x-axis shown 
in Fig.l), I = a/y/2 is the spacing of the bond lattice, 
and Nb is number of sites in the bond lattice (Nt, = 27V). 
The superconducting transition temperature T c (x) is the 
highest temperature at which p s {x, T) is nonzero. We use 



this fact to obtain T c (x) in single-site and cluster mean- 
field theories (the relevant details are summarized in Ap- 
pendix A). As mean field approximations are known [36] 
to overestimate the transition temperature, we treat the 
effect of fluctuations in the model of (Eq.(2)) through MC 
simulations. In these simulations, the standard Metropo- 
lis sampling scheme [51] has been used for planar spins 
{S m = (A m cos 4> m , A m sin(/) m )}, whose lengths are con- 
trolled mainly by J-q (Eq.(2b)). Simulations have been 
carried out for a 100 x 100 square lattice (bond lattice) 
with periodic boundary condition. Typically, 10 5 MC 
steps per spin have been used for equilibration and mea- 
surements were done for next 3 x 10 5 (6 x 10 5 in some 
cases) MC steps per spin. Simulations were done for the 
doping range — 0.4 at various temperatures. 

In our two-dimensional model, true long-range order is 
destroyed by thermal fluctuations, but there is nonzero 
supcrfluid stiffness due to vortex-antivortex binding (the 
BKT transition [39-41]) below a temperature Tbkt- We 
calculate the superfluid stiffness in the MC simulation 
using the formula of Eq.(4) and use it in conjunction 
with the Nelson-Kosterlitz criterion [52] 



Ps(Tbkt) _ 2 
Tbkt t 



(5) 



based on the BKT theory to obtain the vortex binding 
temperature Tbkt{x), which is identical to T c (x) in 2D. 
The above criterion, appropriate for a fixed length XY 
model or cquivalently a low fugacity 2D vortex gas, might 
not give an accurate estimate of Tbkt for the model of 
Eq.(2) in the extreme ovcrdoped regime close to x — x c 
due to large fluctuations in the magnitudes A m [53]. 
Tbkt obtained using Eq.(5) should presumably be quite 
accurate in the underdoped and optimally doped regions 
where the magnitudes effectively become 'frozen' since 
T*(x) » T c (x) resulting in a description of the model 
(Eq.(2)) in terms of an effective fixed-length XY model 
(Appendix C) close to the superconducting transition. 
These results are shown in Fig. 3. Results for the temper- 
ature dependence of the superfluid stiffness are presented 
in Section IV. 

The calculated T c curve is approximately of the same 
parabolic shape as that found experimentally. The causes 
for the qualitative disagreement at both ends (see Fig. 2) 
are not difficult to understand. For very small x, as well 
as for x near x c , our free energy functional needs to be 
extended by including quantum phase fluctuation effects. 
For such values of x, zero point fluctuations are impor- 
tant because the phase stiffness is small. Additionally, 
low-energy mobile electron degrees of freedom need to be 
considered explicitly for x near x c . To include quantum 
phase fluctuation effects, wc supplement the GL func- 
tional of Eq.(2) with the following term that describes 
quantum fluctuations of phases (<p m ) at a minimal level 
[54-56]: 



^Q({Qm}) = 7? ^ QmVmnQn 



(6) 
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FIG. 3: Doping dependence of different temperature scales 
(T° and Tbkt) and the zero temperature gap Ao (Eq.(8b)) 
are shown in the main plot. Inset: Comparison of the T c 's 
obtained from single-site mean-field theory and cluster mean- 
field theory (T c mf and T c cmf respectively) (see Appendix A) 
with the BKT transition temperature Tbkt obtained from 
MC simulation, as discussed in the text. 




FIG. 4: (a) Effect of quantum fluctuation on T c {x) curve of 
Fig.3 for Vb = 0.09Tb. The quantum fluctuation renormal- 
izes T c to T,? throughout the whole x range (Inset). In the 
main figure, we have taken / = 1.55 to change the tempera- 
ture scale To (= 460 K) while keeping A (a; = 0) = 82 meV 
(Section II B) so that the optimal value of Tj* matches that 
of Tbkt in Fig.3. (b) A reasonably good comparison can 
be obtained with experimental T c (x) curve for La214 with 
following choice of parameters (Section II B): x c = 0.345, 
c = 0.33, b = 0.155, / = 1.063, T p = To and Vb = 0.15T with 
A (s = 0) = 82 meV. This choice implies T = 400 K. The 
dip of the experimental T c around x ~ 0.12 is due to the 1/8 
'stripe anomaly' [58] which is out of the scope of the present 
functional of Eq.(2) (see discussion in Section IX). 



Here q m is the Cooper pair number operator at site m, 
and 4> m in Eq.(2c) should be treated as a quantum me- 
chanical operator <j> m , canonically conjugate to q m so that 
[<7m, 4>n] = iSmn [57]. We take the simplest possible form 
for V mn i.e. V mn = Vo6 mn for the purpose of demon- 
strating the effect of quantum fluctuations on the T c (x) 
curve (Fig. 4), where Vq is the strength of on-site Cooper 
pair interaction. We have obtained a single-site mean 
field estimate of T c (x), namely T^(x), including the ef- 
fect of Tq as shown in Fig. 4 and discussed in Appendix 
A. As it is well known, mean field theory overestimates 
the value of the transition temperature. Hence to com- 
pare T®(x) with Tbkt 0*0 of Fig.3 as well as with the 
experimental T c (x) curve, we scale the T® calculated us- 
ing Eq.(A6) by a factor ~ 0.6 in Fig. 4. This factor has 
been estimated by calculating the ratio Tbkt(x)/T™ { (x) 
from Fig.3 (inset). Quantitative agreement for T c for a 
specific cuprate, La2_ x Sr x Cu04 is possible with a partic- 
ular choice of parameters as shown in Fig. 4(b). In this 
extension of the model, we have ignored the long-range 
nature of the Coulomb (or charge) interactions, as well 
as Ohmic dissipation. It has been argued [55] that these 
two factors together result in a fluctuation spectrum sim- 
ilar to the one obtained in an approximation that ignores 
both, but retains the short-range part of the charge in- 
teraction. 

In the remaining parts of the paper, we do not con- 
sider quantum phase fluctuations since they modify the 
results qualitatively only in the extremely underdoped 
and overdoped regions by aborting the superconducting 
transition as the phase stiffness p s (0) becomes small (see 
Fig. 5) at these two extremes in our model. In the rest 
of the x range, these effects are expected to renormal- 
izc [59] the values of the parameters of the functional of 
Eq.(3). We assume that such renormalizations are im- 
plicit in our choice of the parameters A, B and C in tune 
with experimental facts (see Section II B). 



IV. SUPERFLUID DENSITY p 3 {x,T) 

As mentioned above, we have evaluated the superfluid 
density p s at finite temperatures using Eq.(4) by MC sim- 
ulation of our model (Eq.(2)). The results are discussed 
below along with mean-field results. As we have men- 
tioned in Section III, the transition temperature Tbkt 
can be estimated from the universal Nelson-Kosterlitz 
jump of Eq.(5), where p s (T) = above T c . We show 
the results for finite temperature superfluid density in 
Fig.5(a). 

The zero temperature superfluid density can be calcu- 
lated easily from the ground state energy change due to 
a phase twist (a 'spin wave') and is given by 



Ps (x,0) = CA 2 (x) 



(7) 



where Aq(x) is obtained from Eq.(8b) (see Section V). 
Evidently, p s (x, 0) ex x for small x (as is implicit in the 
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timal p s {x, 0) appear, in general, at two different values 
of doping (x C2 /2 > x opt for the present choice of param- 
eters). A similar behavior is observed in experimental 
studies of muon-spin depolarization rate, a a oc p s (x,0) 
of some cuprates which can be sufficiently overdoped 
[17, 18]. The depolarization rate depends on the local 
magnetic field at the location of the muon; this has been 
shown to be proportional to the supcrfiuid stiffness which 
controls the magnetic response of the supcrfiuid [60] . We 
also plot T c (x) as a function of p s (x,0) ('Uemura plot', 
inset of Fig.5(b))which compares well with experimental 
plots of T c vs. (To, measured at low temperatures and 
shown in Refs.17, 19. 

At low temperatures the calculated p s {x,T) decreases 
linearly with T from its zero temperature value i.e. 
p s (x, T) = p s (x, 0) — p' s (x) T; the coefficient of the linear 
term, namely p' s (x) remains more or less independent of 
x for small x and approaches a constant value as x — > 
on the undcrdoped side. The same trend can be observed 
in the experimental data [20, 21] for in-plane magnetic 
penetration depth A Q b, where A~ b oc p s . It is interesting 
that a model for superconductivity such as ours, which 
does not explicitly include electron degrees of freedom 
leads to a linear decrease [61, 62], in the light of the fact 
that the linear dependence has been attributed to ther- 
mal, nodal quasiparticles of the d-wave superconductor 



FIG. 5: (a) Calculated finite temperature superfluid density 
for different x values. The dashed line corresponds to the 
size of universal Nelson-Kosterlitz jump (Eq.(5)) expected 
at a BKT transition. Tbkt(:e) has been obtained from the 
intersection of this line with p 3 (x,T) vs. T curves. Inset: 
p's(x), estimated by fitting p„(x,T) vs. T with a linear form, 
ps(x,T) = p 3 (x,0) — p' s (x)T. (b) Zero temperature superfluid 
density p s (x,0), as a function of x, compared with Tbkt(^) 
and T'c mf ( x )- The superfluid density has been expressed in 
units of energy (meV) as appropriate in 2D. Vertical dashed 
lines indicate x's corresponding to optimal values of p s (x,0) 
and Tbkt(^). The inset shows the 'Uemura plot' [16, 17, 19], 
T c (x) vs. p 3 (x,0). The initial part of the upper branch corre- 
sponds the underdoped region, where the Uemura relation was 
inferred [16] originally. The subsequent decrease of p s (x,0) 
along with T c in the overdoped regime (lower branch) is ob- 
served for example in T^I^CuOe+a [17, 19]. 



choice of C). T c (x), of course, is also proportional to x 
for small x, as can be easily verified from Eq.(A5) (see 
Appendix A), which gives a quite accurate estimate of 
T c for low hole doping. Hence, the Uemura relation [16] 
is seen explicitly to be satisfied for this choice of C. In 
Fig. 5(b) we plot p s {x, 0) as a function of x along with 
T c {x). p s {x, 0) initially increases with x to reach a maxi- 
mum value (so that p y ' = 0) slightly on the overdoped 
side at x = x c . 2 /2 and then ultimately drops to zero at x C2 
as T c also does (see Fig. 3), but the optimal T c {x) and op- 



AVERAGE LOCAL GAP A(x,T) AND THE 
PSEUDOGAP 



The energy gap A m is a thermodynamic variable with 
a certain probability distribution given by the functional 
of Eq.(2). There is no direct measurement of the energy 
gap, unlike that of T c or of the supcrfiuid stiffness dis- 
cussed in Sections III and IV. The information about 
the energy gap is obtained via the coupling of the gap 
(or more precisely, of electron pairs giving rise to the 
gap) to electrons, photons, neutrons etc. In this sec- 
tion, we compute the thermodynamically averaged local 
gap A(x,T) =< A m > and compare our results with 
the broadly observed trends for gaps as inferred from a 
number of measurements on a variety of cuprates. These 
trends are for the pseudogap as a function of hole doping 
x, and for the ratio of the zero temperature gap to the 
pseudogap temperature T*(x) as well as to the directly 
measured superconducting T c . 

Fig. 6 shows the dependence of A(cc,T), calculated in 
single site mean field theory (see Appendix A), on tem- 
perature for different values of the hole doping x. We 
have checked that the values of (A m ) obtained from MC 
simulations are quite similar to the mean-field results, the 
main difference being that the singularity of the mean- 
field values at T c (x) is smoothed out in the MC results. 
Note that the quantity A m = \ip m \ is not the order pa- 
rameter for superconductivity and its average A(cc, T) 
can be (and is) nonzero at temperatures above T c . The 
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FIG. 6: Panel (a) shows the onset of second gap feature in 
A = (A m ) at T c due to the presence of the C term in Eq.(2). 
The dashed lines compares A = (A m )o with A(see text). 
Panel (b) compares the temperature dependence of A for 
T p = To (solid lines) and for T p = 0.65T (dashed lines). 
A changes much more rapidly, especially in the underdoped 
side, with decreasing temperature across T^{x) for the sec- 
ond case. The results shown here and in Fig. 7 were obtained 
from single-site mean-field theory. 
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FIG. 7: (a) Extraction of T ms (x) from the positions of the 
maximum of | g™ | = |f9r(A m )| vs. T curves (upper panel) 
at various doping values. Two local maxima appear in the un- 
derdoped regime, one sharp peak at T c and a broad maximum 
at Tms- Tm S (x) merges with T c (x) in the overdoped side (in- 
set of upper panel) . Similar analysis (lower panel) is carried 
out on | 8 y | (see text for definition) to extract T^ s . (b) 
Comparison of T*(x), identified with T^ 1 , with other rele- 
vant temperature scales; different pseudogap scenarios [7] are 
naturally embodied in our results, as discussed in the text. 



average gap increases smoothly as T decreases; the in- 
crease can be rather abrupt or gradual, depending on the 
parameters (see Fig. 6(b)). The part in A(x, T) 'turning 
on' at T c is generally small. The zero temperature gap 
A (a;) = A(x, 0), is the sum of these two, a gap which 
would have been there even in the absence of phase co- 
herence (shown by the dotted line and calculated from 
A = (A m )o, where the thermal average is evaluated us- 
ing the single site term To of Eq.(2)) and another, due 
entirely to phase coherence. 

Measurements detect a diminution in the density of 
electron states, one which depends on the direction of 
k along the Fermi surface. Different measurements 
(e.g. NMR, resistivity, ARPES etc.) show character- 
istic changes at temperatures which differs by 20 K to 40 
K [5]. The 'pseudogap temperature' T*(x) is, therefore, 



not very well-defined. T* is generally seen to decrease 
with hole doping x, nearly linearly, till it 'hits' the T c (x) 
curve, around (but slightly beyond) x opt . What happens 
next is a matter of considerable controversy. Broadly, 
three scenarios have been argued for, as described for 
example in Ref.7. One of them [63] suggests that the 
pseudogap temperature merges with T c (x) a little beyond 
optimum doping. Another scenario [8, 37, 38] is that it 
goes through the T c (x) dome, reaches zero at a putative 
quantum critical point x qcp , which controls the univer- 
sal low temperature behaviour of the cuprate around it 
in the (x,T) plane. A third [7] is that there is no T* 
beyond the hole concentration x\ at which it 'touches' 
T c (x). Operationally, we identify the pseudogap temper- 
ature as one at which the absolute value of the slope of 
A(x, T) as a function of temperature is a local maximum, 
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calling it T ms (x). In general, this definition leads to two 
characteristic temperatures. One of them is at T c because 
a part of A(x, T) suddenly turns on at T c due to the on- 
set of global phase coherence, leading to a divergence of 
the temperature derivative at T c . The other is at a tem- 
perature higher than T c (x) till an x value slightly above 
x op t- This fact leads to two kinds of behaviour for T ms (x) 
(Fig. 7) and thus for the pscudogap temperature T*{x) if 
these two are identified with each other. If we start from 
the low doping (small a;) side, where T ms (x) is high and 
follow it as x increases, noticing its origin in local pair- 
ing and existence even when there is no global order, 
we see that this branch of T ms (x) denoted as T^ s (x) in 
Fig. 7 hits the T c (x) line at x\ (Fig. 7(b)), goes through 
the T c dome to zero temperature at 'x qcp ' and continues 
to be zero thereafter. On the other hand, if beyond x\ 
we choose the other solution for T ms (x) (called T^ ls (x) in 
Fig. 7), which exists because of the long range order caus- 
ing 'Josephson' or C term in Eq.(2c), then one has a pseu- 
dogap curve which is above T c (x) till x\ and is the same 
as T c (x) thereafter. These are two of the pseudogap cat- 
egories mentioned above. Different types of experiments 
are likely to probe different types of pscudogap. For ex- 
ample, if superconducting phase coherence is destroyed 
with a magnetic field, so that the C or Josephson term 
is ineffective, the observed pscudogap behaviour with x 
is that of the first category. 

At zero temperature the phase coherent classical 
ground state can be represented in terms of nearest- 
neighbor singlet bond pair fields ip m or equivalently i/>i M 
(see Fig.l) as 
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Here, A (x) is the zero temperature gap (see Fig. 3), 
Ao(0) = l/(fVb) and x C2 = x c /(l — 2cx c ) is obtained 
from A(x C2 ,0) - 2C(x C2 ) = 0. 

Our choice of the values of b and / fixes the ratio 
2Ao/To = 2/(/\/&) to be around 3—5, which implies 
that 2Ao(x)/T*(x) also stays close to these values in the 
underdoped regime (Fig. 8). It has been widely reported 
[11, 24] that the ratio of the low temperature ('zero 
temperature') gap to the pseudogap temperature scale, 
specifically A (x)/T*(x), for a range of hole doping, espe- 
cially below the optimum x, is about 4.3/2, which is the 
universal d-wave BCS value [64] for the ratio of zero tem- 
perature gap to superconducting transition temperature. 
Further by choosing c = 0.3, the ratio 2Ao(x)/T c (x) near 
optimal doping is see to be around 10 to 15, as observed 
in cuprates [9, 11], being substantially higher than the 
BCS ratio. In Fig. 8, the ratio 2Ao(x)/T c (x) is shown to 
be more or less constant around optimal doping. The 
increase of this ratio as (1 — x/x C2 )~ x / 2 for large x is an 
artifact of the chosen classical functional. 




FIG. 8: 2A (x)/T c (x) and 2A (x)/T*(x) as functions of x. 
Here T*(x) refers to T^ s (x) (see Fig.7). The long-dashed line 
corresponds to the nearly constant value of 2Aq(x)/T c (x) near 
optimal doping. 



VI. SPECIFIC HEAT 

The electronic specific heat of the superconducting 
cuprates has been measured in many experiments [26- 
28] . It consists of a sharp peak near the superconducting 
transition temperature T c (x) and a broad hump around 
the pseudogap T*(x) [29], both riding on a component 
that is clearly linear in T at temperatures T > T* in op- 
timally doped and overdoped samples. Here, we summa- 
rize theoretical results for the specific heat arising from 
our functional (Eq.(2)), both with and without magnetic 
field. A detailed description is given in a separate paper 
[25]. The functional captures the thermodynamic prob- 
ability of (bosonic) Cooper pair fluctuations and yields 
the contribution of these fluctuations to the specific heat. 
Because of our use of a classical functional, the low tem- 
perature behaviour dominated by quantum effects is not 
properly accounted for; we discuss this below. The low 
energy electronic degree of freedom ignored in our treat- 
ment are the fcrmionic, non-Cooper-pair ones of the de- 
generate electron gas. We use the free energy functional 
(Eq.(2)) to write the specific heat as 
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where §£ = {f 2 exp(T/T p ) + A/T p ) for the particular 
choice of A as in Eq.(3a). Clearly the second term in 
Eq.(9) arises from the fact that T is an effective low 
energy functional whose basic parameters, e.g. A, can be 
temperature dependent. We evaluate C v from Eq.(9) for 
different values of doping x and temperature T by MC 
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FIG. 9: (a) Specific heat obtained from MC simulation of our 
model (Eq.(2)). Panel (b) shows the evolution of the broad 
maximum around T* with doping in the underdoped region. 



sampling of finite 2D systems as mentioned in Section 
IV. The simulations have been carried out with / = 2 
(see Section II B) while choosing A (a; = 0) ~ 54 meV, 
so that T = 400 K. 

We notice that in both theory (see Fig. 9) and 
experiment [28, 65, 66], there is a sharp peak in C v around 
T c (or Tbkt in our case to be more precise). The peak 
amplitude increases as x increases, leading to a BCS like 
shape in the overdoped side. In addition, there is a hump 
[29], relatively broad in temperature, centered around 
T* . The hump is most clearly visible in the calculation 
for the underdoped regime where T* and T c are well sepa- 
rated; its size in the theory depends on A and B (Eq.(3)). 
In experiments, for the underdoped side, its beginnings 
can be seen; unfortunately there are very few experiments 
over a wide enough temperature range to encompass the 
hump fully in this doping regime. The two features, 
namely the peak and the hump, and their evolution with 
x can be rationalized physically. The peak is due the low- 
energy pairing degrees of freedom which cause long-range 
phase coherence leading to superconductivity; these are 
phase fluctuations in the underdoped regime. The hump 
is mainly associated with the regime where the energy 
associated with order parameter magnitude fluctuations 
changes rapidly with temperature. Since this change is a 
crossover centered around T* rather than a phase transi- 
tion, there is only a specific heat hump, not a sharp peak 
or discontinuity. For small x, T* » T c and so we see 
that the hump is well-separated from the peak. As x in- 
creases, T* approaches T c , and in the overdoped regime, 
these are not separated, and there is no hump, only a 
peak corresponding to the superconducting transition. 

In order to compare our results with experiments, in 
particular the features related to critical fluctuations near 
T c ), we remove the contributions that are special to the 
chosen classical functional and are not connected with 
the Cooper-pair degrees of freedom in the real systems. 



FIG. 10: (a) The 'critical' peak appearing near T c for three 
values of x. The inset demonstrates the procedure used for 
the subtraction of the 'non-critical' background (dashed line), 
as mentioned in the text, (b) Analogous plot for the ex- 
perimental specific heat data for Yo.8Cao.2Ba2Cu307_a from 
[27]. Here, x values are estimated using the empirical form 
of Persland et al. [67] . Again, the inset shows the subtracted 
background (dashed line) for x = 0.15. 



Firstly, at low temperatures, T « T c , the fact that we 
have a classical functional here leads to a large specific 
heat of the order of the Dulong-Petit value and there 
is an additional contribution (ex ^, see Eq.(9)) due to 
temperature dependence of A, whereas the actual spe- 
cific heat is expected to be small because of quantum 
effects (it is ~ T 2 due to nodal quasiparticlcs [68]). To 
account for this difference, we compute the leading low- 
tempcrature contribution to the specific heat arising from 
our functional (Eq.(2)). Similarly at high temperatures 
T > T*, the contribution from pairing degrees of freedom 
for the actual system is expected to be small, whereas 
from the functional (Eq.(2)) it is not so due to the sim- 
plified from used for the single-site term (Eq.(2b)). We 
compute C v from a high temperature expansion for the 
intersite term in Eq.(2). We interpolate for the specific 
heat using the low and high temperature expansion re- 
sults, and subtract the resulting part (includes the hump) 
from the calculated specific heat. This subtracted spe- 
cific heat is plotted in Fig. 10(a) for three values of doping. 
These are compared with the experimental electronic spe- 
cific heat data of Rcf. [27] for YBCO after analogous 
subtraction of a 'non-critical' smooth part obtained from 
interpolation between low and high temperature regions 
(excluding the peak) is done (see inset of Fig. 10(b)). This 
procedure also removes linear T contribution to specific 
heat arising from unpaired low energy electronic degrees 
of freedom present in the system but not in our functional 
(Eq.(2)). Since the peaks are large and occur over a nar- 
row temperature near T c , they are relatively free from 
possible errors due to the subtraction procedure men- 
tioned above. The experimental and theoretical results 
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visible for small x, and occurs even for magnetic fields 
as small as a few Tesla. We assume that only the in- 
tcrsite term depends on the vector potential A, via the 
Peierls phase factor, namely that (<p m — </>„) in Eq.(2c) 
is replaced by (<f) m — <j) n — || JL A.dl). The resulting 
specific heat 'peak' curves obtained from MC simulations 
are plotted in Fig. 12 for two x values at different values of 
fn = H1 2 /Qq i.e. the flux going through each elementary 
plaquette of the bond lattice in units of the fundamen- 
tal flux quantum $o = hc/(2e), where H is the applied 
uniform magnetic field perpendicular to the plane (i.e. 
H = Hz) and we assume the extreme type-II limit. The 
results compare well with those of experiment [31]. 



FIG. 11: (a) Evolution of the height of the specific heat 
peak appearing near T c with doping, compared with the 
analogous plot (b) obtained from experimental data for 
Yo.8Cao.2Ba 2 Cu30 7 _ 5 [27]. 
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FIG. 12: Effect of a magnetic field on the specific heat peak 
for (a) x=0.11 and (b) x=0.16. The subtraction procedure 
employed in Fig. 10 is used here as well, as shown in the insets. 



for specific heat peaks are shown separately in Fig. 10. We 
see that they compare well with each other. The qualita- 
tive agreement is brought out clearly in Fig. 11 where we 
plot the specific heat peak height with x and compare the 
dependence with what is observed in experiment. This 
implies that our model for the bond pairs and their inter- 
action to generate a d-wave superconductor is a faithful 
representation of the relevant superconductivity related 
degrees of freedom. 

The effects of a magnetic field on the specific heat 
have been cataloged in [30, 31] where it is found that 
the specific heat peak near T c is increasingly smoothed 
out with magnetic field, but the peak position does not 
shift by much, especially in highly anisotropic systems 
such as Bi2212 and Bi2201. This effect is most clearly 



VII. VORTEX STRUCTURE AND 
ENERGETICS 

We use the functional (Eq.(2)) to find the properties of 
vortices that are topological defects in the ordered phase. 
This has been extensively done in the GL theory for con- 
ventional superconductors [32]. We use the free energy 
functional of Eq.(2) at T = 0, where it describes the 
ground state properties, to generate a single vortex con- 
figuration by minimizing T with respect to A m and <f> m 
at each site while keeping the topological constraint of 
total 27T winding of the phase variables at the boundary 
of a Nb x Nb lattice. This is a standard way of generating 
a stable single k = 1 vortex configuration with the vor- 
tex core at the middle of the central square plaquette in 
the computational lattice. The results for {A m ,</> m } are 
shown in Fig. 13 for two different values of hole doping x, 
namely x = 0.10 (underdoping) and x = 0.30 (overdop- 
ing). Fig. 13(a) shows the order parameter at a point m 
on the square lattice as an arrow whose length is propor- 
tional to the value of A m there, and whose inclination to 
the x-axis is equal to the phase angle <j> m . 

We notice that for the underdoped cupratc (e.g. x = 
0.10) unlike the overdoped one (x = 0.30), the order pa- 
rameter magnitude does not decrease by much as one 
moves radially inwards from far to the core (Fig. 13(b)). 
This is characteristic of a phase or Joscphson vortex 
whose properties have been investigated for coupled 
Josephson junction lattice system [69]. We propose there- 
fore that vortices in cuprates in the underdoped regime 
are essentially Josephson vortices. This is natural here 
because the Cooper pair amplitude A m has sizable fluc- 
tuations only close to T* which is well separated from 
T c (T c << T*) in the underdoped regime so that near 
T = 0, there are very small A fluctuations. Further, for 
a lattice system (and not for a strict continuum) such a 
defect is topologically stable since the smallest possible 
perimeter is the elementary square. On the other hand, 
beyond optimum doping where, according to Fig. 7, T* 
coincides with T c , the order parameter magnitude A m 
decreases substantially on moving radially inwards to- 
wards the vortex core, very much like a 'conventional' 
superconducting or BCS vortex. The variation of the 
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FIG. 14: (a) The excess energy of a vortex AE V as a function 
of system size (see main text) for three values of x. Intercepts 
of the dashed lines with the vertical axis yield the values of 
the corresponding core energies E c . (b) E c is compared with 
T c . Like p s (0) (see Fig.5(b)), E c peaks at x ~ 0.19. The inset 
shows the proportionality of E c and Tbkt in the underdoped 
side. 



FIG. 13: (a) Single vortex configuration for x = 0.10 and 
x = 0.30. Arrows indicate the equivalent planar spins. A 
sublattice transformation has been performed on the phases 
for convenience of representation, (b) Variation of the magni- 
tude of the bond pair field near the vortex core for the afore- 
mentioned values of x. The magnitude is plotted in units of 
its maximum value attained in the bulk, Ao (mentioned at 
the top of each color bar), (c) The angular averaged gap 
magnitude A(r) (normalized by Ao) as function of distance 
from the core for the two x values. Inset shows the doping 
dependence of the magnitude at the core, A CO ro, estimated by 
fitting A(r) with Ao tanh (r/£ c ) + A CO re, while £ c and A core 
are kept as fitting parameters. 



increasing hole density x. 

The core energy E c of a single vortex is naturally de- 
scribed as the extra energy AE V = E v — Eq where Eq is 
the energy of the ground state configuration (the Ne'el 
ordered state in this case) and E v is the total energy of a 
single vortex configuration, from which the elastic energy 
due to phase deformation [36] is subtracted, i.e. 



A£L 



E c + ir Ps (0)ln{R/l) 



(10) 



normalized magnitude of the bond pair field A(r)/Ao 
with the radial distance r from the vortex core in the 
two cases is shown in detail in Fig. 13(c), which clearly 
illustrates the difference between the behavior in the two 
cases. The inset of Fig. 13(c) shows the extrapolated val- 
ues of the magnitudes (A core ) at the core (r = 0) as a 
function of x, indicating that there is a smooth crossover 
from a Josephson-like vortex to a BCS-like vortex with 



The quantity R is defined as R = (Nf,—l)l/-\/%, where I is 
the lattice constant of the bond lattice, so that ttR 2 is the 
area of the computational lattice. We plot in Fig. 14(b) 
the core energy E c as a function of x, both its absolute 
value and its ratio with T c . E c has been estimated from 
the intercept of the AE V vs. hx(R/l) (different system 
sizes) straight line with the energy axis. We notice that 
for small x, E c (x) ex T c (x) (inset of Fig.l4)(b), not sur- 
prising from XY model considerations [70]. 
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VIII. ELECTRON SPECTRAL FUNCTION AND 
ARPES 

The cuprate superconductor obviously has both elec- 
trons, and Cooper pairs of the same electrons, coexisting 
with each other. In a GL-like approach such as ours, only 
the latter are explicit, while the former are 'integrated 
out'. However, effects connected with the pair degrees of 
freedom are explored experimentally via their coupling 
to electrons, a very prominent example being photocmis- 
sion in which the momentum and energy spectrum of 
electrons ejected from the metal by photons of known 
energy and momentum is investigated. Since ARPES 
(angle resolved photoemission spectroscopy) [9, 10] is a 
major and increasingly high-resolution [71] source of in- 
formation from which the behaviour of pair degrees of 
freedom is inferred, we mention here some experimental 
consequences of a theory of the coupling between elec- 
trons and the complex bond pair amplitude ip m . The 
theory as well as a number of its predictions (in agree- 
ment with ARPES measurements) are described in detail 
in Ref.35. 

In formulating a theory of the above kind, one faces 
the difficulty of having to develop a description of elec- 
trons in a presumably strongly correlated system such 
as a cuprate, which is viewed as a doped Mott insula- 
tor [1] with strong low-energy antiferromagnetic correla- 
tion between electrons at nearest neighbor sites [4]. In 
particular, one needs to commit oneself to some model 
for electron dynamics which then implies an approach 
to the coupling between electronic and pair degrees of 
freedom. We develop what we believe is a minimal the- 
ory, appropriate for low-energy physics. We assume that 
for low energies |w| < Ao, well-defined electronic (tight- 
binding lattice) states with renormalized hopping ampli- 
tudes t, i', t" etc. exist and couple to low-energy pair 
fluctuations ip m = ip ilt = ((a^Oi+^t - a i - t a i+f _ li )/2) (sec 
Fig.l). Superconducting order (more precisely, phase 
stiffness) and fluctuations in it are reflected respectively 
in the average (Vv( T )) ana - the correlation function 
(^{t)^*^(t')) (or its Fourier transform D fJ , fJ ,'(2q,iz m ), 
z m = 2rmr/f3 being the bosonic Matsubara frequency 
where m is an integer). A nonzero value of (ipifi (t)) in the 
'AF' long-range ordered phase below T c leads to the well 
known Gor'kov d-wave Green's function and quasiparti- 
cles with spectral gap Ak = (A^/ '2) (cos k x a — cosk y a). 
The correlation function D^i (q, oS) has a generic form 
for small q and u which can be related to the functional 
(Eq.(2)). 

The coupling between low excitation energy electrons 
and low-lying pair fluctuations (both inevitable) leads to 
a self energy with a significant structure as a function of 
electron momentum k and excitation energy ui. Physi- 
cally, we have electrons (e.g. those with energy near the 
Fermi energy) moving in a medium of pairs which have fi- 
nite range 'AF' or d-wave correlation for T > T c and have 
long-range order of this kind for T < T c (in addition to 
'spin wave' like fluctuations). The electrons exist both 



as constituents of Cooper pairs and as individual enti- 
ties; the pairs and the electrons are in mutual 'chemical' 
equilibrium. The energy shift or dynamic polarization of 
electrons due to this process leads to a number of effects 
which are described in [35]. For example, for T > T c we 
find a pseudogap in electronic density of states which per- 
sists till T*. We get Fermi arcs [9, 10, 72] i.e. regions on 
the putative Fermi surface where the quasiparticle spec- 
tral density has a peak at zero excitation energy in con- 
trast to the pseudogap region where the peak is not at the 
Fermi energy. The antinodal pseudogap 'fills up' between 
T c and T* with increasing temperature. Below T c , there 
is a sharp antinodal quasiparticle peak whose strength 
is related to the superfluid density as observed in exper- 
iment [73]. We also obtain a 'bending' or departure of 
the Ak vs. k curve from the mean-field canonical d-wave 
form due to order parameter or 'spin wave' fluctuations. 
Here we only outline our theoretical approach and show 
how a temperature T an can be obtained from the filling 
in of the antinodal pseudogap above T c . We find that 
T an compares well in its magnitude and x-dependence 
with other measures of the pseudogap temperature scale 
described in Section V. 

The physical quantity of interest is 



A(k,i 



-Im[G(k, iv n ->• u + iS) (11) 



(the fermionic Matsubara frequency, v n = (2n + l)n/j3, 
n being an integer). Assuming translational invariancc 
one has the Dyson equation for G, namely 

G-^ivn) = (G°)- 1 (k,z^)-S(k,z^) (12) 

where E(k, iv n ) is the self energy. 

G (k, iv n ) is described in terms of a spectral density 
in the usual Lehmann representation [74]. The spectral 
density for low excitation energies has a Dirac (5-function 
part i.e. A (\s.,U)) = Zk5(uj — £k) where £k is the ef- 
fective quasiparticle energy measured from the chemical 
potential jj, and Zk (< 1) is the quasiparticle residue. 
In the 'plain vanilla' or renormalized tight-binding free- 



MI 



particle theory [75, 76] Zk = 1 and £k 
with e k = 9t SfRi-R •) *tj : exp[ — ik.(Rj — Rj)], so that 
G°(k, iv n ) — l/{ii> n — £k)- The factor g t is due to corre- 
lation effects calculated in the Gutzwillcr approximation 
[76] which projects out states with doubly occupied sites; 
one further assumes that the renormalized quasiparticles 
propagate coherently. 

We use a standard approximation for E(k, iv n ) which 
is shown diagrammatically in Fig. 15. This describes a 
'phonon' like process neglecting vertex corrections; the 
propagating electron become a Cooper pair (boson) plus 
an electron in the intermediate state; these recombine to 
give a final state electron with the same (k, iv n ). The 
internal propagator in Fig. 15 is the true or full propa- 
gator G. However, in common with general practice, we 
find E and thence G by inserting G° instead of G in the 
former. This is known to be quite accurate [74], e.g. for 
the coupled electron-phonon system. 
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Vv„ n 



-k+2q, -4v n + 



iz„ 



u> Mv„ 



FIG. 15: Self energy approximation used to calculate the elec- 
tron Green function G(k, iu n )- The wavy line denotes the 
pair propagator D Mfl /(2q, iz m ) and the line with an arrow- 
head pointing towards left indicates the full electron Green 
function G(— k + 2q, — iu n + izm.) (see text). The external 
lines (dashed) at two ends of the diagram represent bare (left) 
and true (right) electron propagators. In the static approx- 
imation £> MM /(2q, izm) = (D l j,f 1 ,>(2q)/T 2 )S Zm ,o and the sum- 
mation over the internal bosonic Matsubara frequency in the 
above diagram drops out (see Appendix B). 



In the static approximation valid at high temperatures 
when the pair lifetime t p » l/(ksT) (see Appendix B), 
the general algebraic expression for £(k, iu n ) is 

E(k,w/ n ) = 

^ J2 G°(-k + 2q,-w n ) J D w /(2q)/ /i (k,q)/^(k,q) 



N 



(13) 



where N is the total number of Cu sites on a sin- 
gle CuC>2 plane and /x, // refer to the direction of 
the bond i.e. x or y. The static pair propaga- 
tor is IV (2q) = T 2 IV (2q,0) (see Fig. 15) where 
£>w'(2q) = ER^^'(R)exp(-i2q.R) with D^(R) = 
(tp^(R,)il>*,(0)). Since the XY-like interaction term 
(Eq.(2c)) between nearest-neighbor bond pairs (see 
Fig.l) is antiferromagnetic, 



D XX (R) = £>„„(R) = -D xy (R) = D(K) 



(14) 



Further, the quantity / M (k, q) is a form factor describ- 
ing the coupling between an electron and a bond pair. 
For a tight binding lattice and nearest-neighbor bonds, 
//*(k,q) = cos[(fc M -g M )a]. 

The pair correlator of Eq.(14) can be written in the 
standard way [36], 

£>(R m - R„) = (VvnXC)) + S(Km - R„) (15) 

where ifj m = A m exp (iip m ) with ip m — <p m for x-bonds 
and ip m = cj) m + it for y-bonds (sec Fig.l); 5(R) is 
the fluctuation term. In the long-range ordered state 
below T c , the first term is nonzero. In that case, if 
one neglects effects of fluctuations i.e. 5(R) altogether 
(as is done in mean-field theory), then one obtains 
the exact Gor'kov self energy form [74] i.e S(k, iv n ) — 
^k/(*^« +60 m Eq.(13) and corresponding spectral gap 



Ak = (Arf/2)(cos k x a — cos k y a) in the Necl ordered state. 
Spin-wave-like fluctuations below T c can be incorporated 
through S(R) which generally decays algebraically for 
large distances i.e. S(R) ~ R~ v (17 > 0, its value de- 
pends on dimension). Above T c , (ip(R)) = and the only 
contribution comes from the fluctuation part. Gencri- 
cally, there is a finite correlation length £ above T c and 
5(R) ~ exp(-R/0 or 5(q) ~ 1/[1 + (£g) 2 ]. 

Since we are mainly interested in the spectroscopic fea- 
tures of the pseudogap regime when T* (x) is perceptibly 
higher than T c (x) so that fluctuations in the pair magni- 
tude A m are small and short ranged, we write, 



D(R) ~ < A(R) >< A(0) > ( e J [v( R )-^(°)]) 
= A 2 F{R) 



(16) 



where F(R) = (e'^-^ )!) is the phase correlator. 

Analytical expression for the self-energy from Eq.(13) 
can be obtained below T c , where quasi- long-range or- 
der in purely 2D system or true long-range order in 
anisotropic 3D system occurs, as well as above T c in the 
temperature regime where the exponential decay of cor- 
relation is governed by a large correlation length £ [35]. 
We have carried out calculations [35] for both anisotropic 
3D and 2D cases, while incorporating a small interlayer 
coupling Cj_ (with C/C± - 100 as suitable for Bi2212) 
in Eq.(2) for the former. Above T c the anisotropic 3D 
system behaves effectively as 2D [77] and our results for 
various spectral properties are quantitatively similar and 
even below T c , for this large anisotropy ratio, qualita- 
tive features are the same for both the cases. Hence, we 
present here the results for the pure 2D system. More 
specifically, here we have used the form 



F(R) = (Kr) e- R ^ 



(17) 



to calculate the self energy (Eq.(13)). Here A is re- 
lated to the upper wave-vector cutoff of the lattice and 
r) = T/(2"Kps) below T c where £ — > 00. Above T c , we 
have set 77 = j^bkt = 0.25. A combination of MC simula- 
tion and well-known Kostcrlitz-Thouless renormalization 
group relations has been used to estimate £(x, T) from 
the functional (Eq.(2)) (see Appendix C for details). The 
self energy S(k, iv n ) obtained using the form of F(R) in 
Eq.(17) evolves smoothly from below T c (superconduct- 
ing state) to above T c (pseudogap state). 

For k on the Fermi surface [78] in the antinodal re- 
gion, we calculate A(k = k an ,w). Above T c but below a 
certain temperature (denoted as T an ), two peaks appear 
in _4(k an ,u;) at nonzero u, one at us < and another at 
uj > 0, signaling the presence of a pseudogap above T c . 
The antinodal gap (denoted as A an ) can be defined from 
the position of the peak at negative energy (w < 0). This 
quantity has been plotted in Fig. 16(a) as a function of 
temperature for a few values of x. 

The quantity A an goes to zero rather abruptly at T an , 
though the average local gap A is non zero above T an (see 
Fig. 6). The antinodal pseudogap fills in at this tempera- 
ture [35]. In Fig. 16(b), T an is plotted as a function of x. 
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FIG. 16: (a) Variation of antinodal gap A an with tempera- 
ture. Slope discontinuities in A an vs. T curves correspond 
to T c (Tbkt). (b) Pseudogap temperature scale T an ob- 
tained from the antinodal gap filling criterion mentioned in 
the main text. T an (x) is compared with other temperature 
scales, T®(x), T^ a (x) and Tbkt(3:). Here, we have taken the 
nearest-neighbor hopping t = 300 meV and the next-nearest- 
neighbor hopping t' — — i/4 [81]. 



We notice that this temperature is close to various pseu- 
dogap related temperatures e.g. the somewhat arbitrary 
linear T®(x) used in Eq.(2), as well as the temperature 
scale T^ s (x) estimated from the temperature dependence 
of the local gap magnitude. The x-dependence of T an is 
similar to that of T* as inferred from ARPES [79] as well 
from various other probes such as Raman spectroscopy 
[12] and spin susceptibility [5, 8] over a rather large range 
of x. 

The picture used in our calculation continues to regard 
the electrons as coherent at all temperatures whereas 
there is experimental evidence [80] that the incoherence 
temperature is proportional to x so that it is rather small 
for small x. Also, for very small x, the holes tend to lo- 
calize, so that a renormalized band theory implying ex- 
tended homogeneous electronic states is inappropriate. 



IX. DISCUSSION AND FUTURE PROSPECTS 

We mention here some obvious directions in which the 
functional and the approach used here need to be de- 
veloped. One is to obtain other testable/experimentally 
measured consequences of the proposed functional. For 
example in a magnetic field, the intersite term in Eq.(2) 
has its phase altered by Peierls phase factor, as we have 
mentioned at the end of Section VI. One should use this 
to find the T C (H) curve for different values of doping x 
and thence the 'bare' coherence length £o defined through 

the phenomenological equation, -!=- (fjf) T _ T = ( $ )• 

The charge related response of a system described by 
Eq.(2), e.g. the diagonal and off-diagonal components of 
the conductivity tensor, a xx (H 7 T > T c ) and a xy (H,T > 
T c ), and the Nernst coefficient a xy (H,T > T c ), needs to 
be calculated and compared with experiment. Slightly 
farther afield, the coupling of the field ip m to different 
probes will enable one to analyze experimental results 
obtained e.g. from scanning tunneling spectroscopy, Ra- 
man spectroscopy and neutron scattering. A generaliza- 
tion to a quantum ijj m functional and inclusion of other 
time-dependent effects, e.g. Coulomb interaction and 
dissipation may enable one to describe quantum phase- 
fluctuation effects, which are specially prominent (and 
decisive) for extreme underdoping [82]. 

A very peculiar feature of cuprates is the unusually 
large proximity effect [83] observed in them. While XY- 
spin-likc models have been proposed for this [84] , a com- 
plete understanding of the size, temperature and doping 
dependence etc. does not exist. It is possible that the 
present theory can be adapted to address this question. 

The theory presented needs to be extended in many 
major ways. For example, there is a lot of experimental 
evidence [4] that the system is a Mott insulator at x = 0, 
with a large superexchange Jy ~ 0.15 cV., as well as for 
low-energy magnetic correlations in doped cuprates. This 
antifcrromagnetic interaction evolves into superconduc- 
tivity for surprisingly small hole doping, x > 0.05. While 
the crossover and the possibility of coexistence have been 
investigated at T = [85-87] , there is need for a coupled 
functional for these two bosonic degrees of freedom that 
goes over to the kind of theory we have described above at 
large x, while it describes an antiferromagnetic Mott in- 
sulator at x — and persistent spin correlations (includ- 
ing spin density wave correlations) at x ^ 0. Similarly 
there is considerable evidence for other kinds of correla- 
tions, e.g. nematic [88], stripes [89], checkerboard [90], 
and charge density wave [91] whose significance varies 
with material, doping (including commensuration effects 
[58]) and temperature. An appropriate GL like functional 
is one way of exploring the details of this competition: 
attempt in this direction already exist [92]. 

The cuprate properties are very sensitive to certain im- 
purities e.g. Zn replacing Cu. Whether this can be de- 
scribed well in a GL like theory is an interesting question. 
The effect of impurities or in-plane/intra-planc disorder 



17 



is an even more general question in terms of its effect 
on pairing degrees of freedom as well as incorporation of 
this effect in a this kind of picture. A subject of basic 
interest in cuprate superconductivity is the possibility of 
time-reversal symmetry breaking associable with T* [38] . 
There are at least two observations, one of Kerr effect 
[93] and another of ferromagnctism with lattice symme- 
try [94], which seem to point to time reversal symmetry 
breaking below T* . Since these involve spontaneous long- 
range order in circulating electric currents, each within a 
single unit cell of the lattice, and these currents can be 
modeled in a GL functional, one can explore this novel 
phase and its consequences in our theory. 

In conclusion, we believe that the phenomenological 
theory proposed and developed here not only ties to- 
gether a range of cuprate superconductivity phenomena 
qualitatively and confronts them quantitatively with ex- 
periment, but also has the potential to explore meaning- 
fully many other phenomena observed in them. 
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Appendix A: Mean Field Theory 

We describe here various approximate solutions for the 
properties of the lattice functional (Eq.(2)). The approx- 
imations discussed here are single-site mean field theory 
and cluster mean field theory. We also make use of sev- 
eral well-known results from the Berezinskii-Kostcrlitz- 
Thoulcss theory [36, 40, 41] for XY spins in two dimen- 
sions, in combination with Monte Carlo simulation (see 
Section III). For positive C in Eq.('2c), there is a low- 
temperature phase with long range 'AF' order (d-wave 
superconductivity) or broken symmetry (for d > 2). The 
most common approximation for locating and describ- 
ing this transition is (single-site) mean field theory, in 
which we sclf-consistently calculate the staggered 'mag- 
netic field' h = (h x ,h y ), acting on the planar spins 
S m = (A m cos(j) m , A m sin0 m ), due to its nearest neigh- 
bors, assuming it to be the same at each site (modulo the 
sign change due to the two sublattice 'AF' order). 

In such a mean field theory [36, 95], the self-consistent 
solution is given by 



with 



(Sa) 



h a = 4C(S a )o (a = x,y) 



h a \ f™ A*dAP Q (A)h(hA/T) 



h ) f °° AdAP (A)I (hA/T) 



(Al) 



(A2) 



Here, P (A) = exp (-/3(AA 2 + (5/2)A 4 )) dictates the 
local distribution (thermal) of gap magnitude, h = 



h-x + hy is the magnitude of the 'staggered' field and 
Io, I\ are modified Bessel functions of first kind. The 
transition temperature T c (which is denoted as T™ f in 
Fig. 3) satisfies the implicit equation 

2C(A 2 ) Po | T=Tc =T c (A3) 

where (A 2 ) Pf) = J™ A*dAP (A)/ J™ AdAP (A). 
Other physical quantities, such as the supcrfluid stiff- 
ness, the superconducting order parameter, the internal 
energy (and its temperature derivative, the specific heat 
C v ), can be obtained using the self-consistent solution of 
Eq.(Al). For instance, in this approximation, the super- 
fluid density p s is given by 

Ps z _ T^~(y^ A m A m+/i cos(0 m _ (f) m+fl )) 



2N h "^ 



C 2^, w q )oj 



(A4) 



In reality, the field acting on a 'spin' fluctuates from 
site to site. The spatially local fluctuations are system- 
atically included in the well-known cluster theories, the 
oldest of which is the Bcthe-Peierls approximation [95], 
which consists of a single site coupled to the nearest- 
neighbors which are described by a mean field. We have 
used it to calculate an 'improved' T c (T™ f ), as shown in 
the inset of Fig. 3. 

For small x, where amplitude fluctuations can be ne- 
glected, an estimate of T c (denoted as T c> o) is obtained 
by replacing (A 2 )o in the above relation (Eq.(A3)) by 
A 2 n that minimizes the single-site term J-q, so that 
A 2 n = — A(x, T)/B for x < x c and A mj0 = for x > x c . 
In this approximation, 



T, 



i-.O 



2xc / x 

2xc + b \ x c 



x < x r 



= X > X r 



(A5) 



Here we have neglected the exponential temperature de- 
pendence of A (Eq.(3a)). Consequently x op t can also 
be estimated by setting g c x ° = 0, which gives x op t = 
k(VWc) 2 + (2bx c /c)-(b/c)). 

If one includes the term Tq (Eq.(6)), the self- 
consistency condition for T c in Eq.(A3) gets modified in 
the following manner [57], 



UC(A 2 ) P0 J^ 



4C(A 2 ) Po / dT(cos0 m (r)cos(/) m (O))^ Q ] =1 



T=T C 



(A6) 



where the average (...)jr Q is calculated using the eigen- 
states of Tq and the imaginary time on-site phase-phase 
correlator in Eq.(A6) is given by [57] 

1 



rV (l-r//3)_ (A?) 

where Vq is the on-site Cooper pair interaction strength. 



(cos^ m (r)cos^ m (0))jr Q = -e 



Appendix B: Electron Self Energy in Static 
Approximation 

The self energy depicted in Fig. 15 can be written in 
the following form using Go(— k + 2q, —iv n + iz m ) for the 
internal electron propagator, 



S(k, ii 



T 2 

— Y 



D(2q,iz m )V(k,q) 



N z ' IVr, 

q,?n 



e 



k-2q 



(Bl) 



where £>(2q, iz m ) = (1/T) $ dr £ R D(R, r ) e -*2 q .R+« m 7 
is the Fourier transform of the time-dependent propaga- 
tor and T-^k, q) = [cos (k x a — q x a) — cos (k y a — q y a)] 2 . 
If the pairs acquire a finite lifetime r p , the pair 
correlator can be represented in terms of the prod- 
uct of the static propagator (Eq.(14)) and a time- 
dependent part as D(R,t) = £)(R)e _ '/ T p so that 
D(2q, iz m ) = (1/T)(e^/^ - l)D(2q)/(iz m + i/r p ). This 
form indicates that pair correlations decay temporally 
with a lifetime r p (one can instead take an oscillatory 
form i.e. D(R,t) ~ cos(£/t p ) but this does not change 
our main conclusion). One can perform the summation 
over the bosonic Matsubara frequencies (z m ) in Eq.(Bl) 
with the aforementioned form of -D(2q, iz m ) and obtain 



2j(k, iz m ) 
1 
T 



1_ ^ £>(2q)P (k,q)((l - eW^)/(£ k _ 2q ) + e^) 



i(i/ n + l/r p ) + Ck- 



2q 
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Here f(ui) = l/(e' 9 " + 1) is the Fermi function. When 
T >> (1/t p ) (also tVi >> (l/r p ) since z/ n ex T) i.e. in- 
verse pair lifetime is much smaller than T, the self energy 
given above would effectively reduce to the form given in 
Eq.(13). 



Appendix C: Estimation of Correlation Length £ 



We estimate r\ = T/(2Trp s ) (below T c ) and £ (above 
T c ) that appear in Eq.(17). As already discussed, we cal- 
culated p s below T c from our functional in Section IV by 
performing MC simulation. Correlation length £ can be 
estimated by fitting obtained p s {x,T) below T c with the 
BKT form, p 8 (x,T) = p B [T-(a:)][l + b(x)^T c (x) - T\ 
with p s (T~)/T- = 2/tt, and b(x) and T c (x) as fit- 
ting parameters. BKT RG relates [96] b(x) to the 
temperature-dependence of £ above T c through £(x, T) ~ 

where bb' = it/2 and Qq is a 



V{x)/y/T-T c {x) 



aoexp 

microscopic length scale of the order of the lattice spac- 
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